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Abstract 

We study the extension of a phenomenologically successful quasiparticle model that describes 
lattice results of the equation of state of the deconfined phase of QCD for T c < T < 4T C , to finite 
quark chemical potential [i. The phase boundary line T c (/-i), the pressure difference Ap(T, /i) = 
(p(T, fj.) —p(T,fj, = 0))/T 4 and the quark number density n q (T, jj,)/T 3 are calculated and compared 
to recent lattice results. Good agreement is found up to quark chemical potentials of order fj, ~ T c . 
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1 Introduction 



The phase structure of QCD at high temperature and non- vanishing baryon chemical potential has been 
subject of intense research in recent years. Heavy-ion collisions at high energies have been and are being 
explored at SPS/CERN and RHIC/BNL Q] in search for signals of the Quark-Gluon Plasma (QGP). 
Large-scale lattice QCD computations at finite temperature have been performed 03 0] , and first 
extensions to non-zero baryon chemical potential appear now to be feasible. It has proven possible to 
trace out the phase boundary line T c (fi) separating the hadronic phase from the QGP phase for Nf = 4 
[HE! Nj — 2 |7| and Nf = 3 [5] flavors of quarks up to quark chemical potentials /x of order T c . First 
numerical results for the QCD Equation of State (EoS), i.e. the pressure p(T, /x) and the quark number 
density n q (T, /i) are also available for Nf =2 + 10 and Nf = 2 JU] . As well as being of intrinsic 
theoretical interest, such studies provide conceptual guidance for current heavy ion collision experiments 
at SPS and RHIC, where the chemical freeze-out occurs at \if. . — 100 MeV, (baryon chemical potential 
Hb ^ 300 MeV) [IT] and n f . . ~ 15 MeV, (/x B ^ 45 MeV) [X2], respectively. 

Systematic perturbative expansions of the QCD equation of state within the framework of thermal field 
theory show bad convergence even for very large temperatures (several times T c ) far beyond the region 
accessible to present experiments |13| . Various techniques, such as dimensional reduction, screened 
perturbation theory or hard-thermal loop (HTL) perturbation theory show better convergence and good 
agreement with lattice results for T > 3T C Various interpretations of the lattice data have been 

attempted in terms of physical quantities, most prominently as the EoS of a gas of massive quark and 
gluon quasiparticles. Their thermally generated masses are based on perturbative calculations carried 
out in the HTL scheme ^1 IH\ ■ This approach has been extended to non- vanishing quark chemical 
potential and good agreement with finite lattice calculations for Nf = 2 + 1 flavors has been found 
|18j . More recently, the QGP has also been described in terms of a condensate of Z$ Wilson lines [B?] 
and by more refined quasiparticle models based on the HTL-resummed entropy and (NLO) extensions 
thereof |2()| . These models have found support from resummed perturbation theory |21j for temperatures 
T > 3T C . However, they have difficulties explaining the dropping of the thermal gluon screening mass in 
the vicinity of the phase transition. An improved quasiparticle model |22| shows the correct temperature 
dependence of the Debye mass and reproduces lattice thermodynamical quantities such as the pressure, 
the energy density and the entropy density very well. The main new ingredient of this model is a 
phcnomenological parametrization of (de)confinement. 

In the present work, this improved quasiparticle model is extended to finite quark chemical potential 
/i. In section 2, a brief review of the quasiparticle model with confinement is given. The extension of 
the model to finite quark chemical potential /x is discussed in detail in section 3. Numerical results are 
presented in chapter 4. The phase boundary line T c (/x) that separates the hadronic from the QGP-phase 
is discussed and the quasiparticle model result is compared to recent lattice simulations. Results for the 
pressure difference from fj, = and the quark number density for various values of the chemical potential 
/x are also presented and compared to recent lattice simulations. A summary is given in section 5. 

2 Quasiparticle model with confinement 

It is possible to describe the EoS of hot QCD at vanishing quark chemical potential /i to good approxi- 
mation by the EoS of a gas of quasiparticles with thermally generated masses, incorporating confinement 
effectively by a temperature-dependent, reduced number of thermodynamically active degrees of free- 
dom. This method is briefly outlined in this section. For a more detailed discussion the reader is referred 
to ref.j22J. 

At very high temperatures, spectral functions for gluons or quarks of the form 8{E 2 — k 2 — m 2 (T)) with 
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m(T) ~ gT are found in HTL perturbative calculations. Here, E is the particle energy, k the absolute 
value of its momentum, m(T) its thermally generated mass and g the QCD coupling constant. As long as 
the spectral function at lower temperatures resembles qualitatively this asymptotic form, a quasiparticle 
description is expected to be applicable. QCD dynamics is then incorporated in the thermal masses 
of the quark and gluon quasiparticles. These thermal masses are obtained from the self-energies of the 
corresponding particles, evaluated at thermal momenta k ~ T: 



m\ = ml q 



m 2 g = ml g 
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Nf is the number of flavors, N c the number of colors. The effective coupling strength G is specified as 



G(T, M = 0) = 



(Jo 



j T 



(3) 



y/UN c - 2Nf 

Setting g = 9.4, (3 — 0.1, the effective masses as given in equations (QJ and (0), approach the HTL 
result at high temperatures. (A small shift 6 = 10 -6 helps fine-tuning at T ~ T c ). Because of the 
existence of a heat bath background, new partonic excitations, plasmons (longitudinal gluons) and 
plasminos (quark-hole excitations) are also present in the plasma. However, their spectral strengths are 
exponentially suppressed for hard momenta and large temperatures and consequently these states are 
essentially unpopulated |23|. The functional dependence of m g (T) on T is based on the conjecture that 
the phase transition is second order or weakly first order which suggests an almost power-like behavior 
(T - T c f with some critical exponent /3 > 0. It is assumed that the pseudo critical form of the 
effective coupling constant given in equation © also provides the correct approximate expression for the 
effective quark mass. This is supported by a non-perturbative dispersion relation analysis for a thermal 
quark interacting with the gluon condensate |24| . 



Close to T c the picture of a non-interacting gas is not appropriate because the driving force of the 
transition, the confinement process, is not taken into account. Below T c , the relevant degrees of freedom 
are pions and other hadrons. Approaching T c from below, deconfincment sets in and the quarks and 
gluons are liberated, followed by a sudden increase in entropy and energy density. Conversely, when 
approaching the phase transition from above, the decrease in the thermodynamic quantities is not 
primarily caused by increasing masses of the quasiparticles, but by the reduction of the number of 
thermally active degrees of freedom due to the onset of confinement. For example, gluons begin to form 
heavy clusters (glueballs), so that the gluon density gets reduced as T c is approached from above. This 
feature can be incorporated in the quasiparticle picture by modifying the number of effective degrees of 
freedom by a temperature dependent confinement factor C(T): 

C(T,fi = 0) = C ([l + 5 c ]-^Y (4) 

The confinement factor is taken to be universal. The parameters Cq , 5 C and C are fixed by reproducing 
the entropy density that results from lattice QCD thermodynamics. Since the results of lattice calcula- 
tions with dynamical quarks are still dependent on the details of the simulations, Co, 5 C and j3 c should 
be finetuned for different lattice calculations. 



For homogeneous systems of large volume V, the Hclmholtz free energy F is related to the pressure p by 
F(T, V) — —p(T)/V. In the present framework of a gas of quasiparticles, its explicit expression reads 
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f'OC 7,4 rj TV7- f POO 7,4 

PiT) = -^J dkC ^ fB ^W k + ^T,J o dkC(T)f D (El)--B(T). (5) 

v g is the gluon degeneracy factor, E g k = ^Jk 2 + m g (T) is the gluon energy, E k — ^Jk 2 + m 2 {T) the 

quark energy, fB{E k 1 ) — (exp((_E?)/T) — 1) the Bose-Einstein distribution function of gluons and 
fn(E%) = (exp((£|) /T) + the Fermi-Dirac distribution function of quarks. The energy density e 
and the entropy density s take the form 

Nr 

r°° 9 at r 00 

£(T) = ^ i dkk2c ( T )fB(E 9 k )K + £ y o dkk 2 C(T)f D (El)El + B(T). (6) 

and 

i/„ Z" 00 „ ffc 2 + m 2 (T) 27V r ^4 f 00 , - ffc 2 + m2(T) 

The function B{T) is introduced to act as a background field. It is necessary in order to maintain 
thermodynamic consistency: p, e and s = dp/dT have to satisfy the Gibbs-Duhcm relation e + p = Ts — 
Tdp/dT. B{T) basically compensates the additional T-derivatives from the temperature-dependent 
masses in p and thus is not an independent quantity. Since B(T) adds to the energy density of the 
quasiparticles, it can be interpreted as the thermal vacuum energy density. The entropy density, as a 
measure of phase space, is unaffected by B{T). 

3 Finite chemical potential 

The quasiparticle model reviewed in the previous section accurately reproduces lattice thermodynamical 
quantities such as the pressure, the energy density and the entropy density in the temperature range 
T c <T < 4T C at vanishing chemical potential [22] • However, many physical questions, e.g. the structure 
of quark cores in massive neutron stars, the baryon contrast prior to cosmic confinement or the evolu- 
tion of the baryon number in the mid-rapidity region of central heavy-ion collisions, require a detailed 
understanding of the EoS at non-vanishing quark chemical potential. In this section, a thermodynami- 
cally self-consistent extension of the quasiparticle model to finite quark chemical potentials is presented. 
Results for various observables are then computed and compared to finite /i lattice results in the next 
section. 

At vanishing quark chemical potential, it is conjectured from asymptotic freedom that QCD undergoes 
a phase transition from the hadronic phase to the QGP phase. At extremely high density, cold quark 
matter is necessarily in the Color-Flavor-Locked (CFL) phase in which quarks of all three colors and all 
three flavors form cooper pairs. It is expected that this phase is separated from the hadronic phase by 
the color superconducting 2SC phase. For a review of the QCD phase diagram, the reader is referred to 
|25| . Our extension of the quasiparticle model provides a straightforward way to map the EoS at finite 
temperature and vanishing quark chemical potential into the T — fi plane without further assumptions. 
However, since this continuous mapping relies on quark and gluon quasiparticles, it cannot provide infor- 
mation about other possible phases with a different (quasiparticle) structure. It is therefore applicable 
in a limited range of not too large chemical potentials. 

The pressure of an ideal gas of quark and gluon quasiparticles with effective masses depending on tem- 
perature and quark chemical potential, is given by 
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p( T> ri = -JL dAC f(T, A i)/ B (Eg)- ? + -|£/ &C(T, / .)[/+(^)+/ D (^)]- ? -B(T l/1 ), (8) 

with f^{El) = (exp((i?*=F/i)/T) + l) _1 . The effective coupling strength G(T, /j), the confinement factor 
C(T,ii) and the mean field contribution B(T,[i) now also depend on the quark chemical potential /i. 
B(T,ii) is calculated in appendix 1X1 The quark number density (which is related to the baryon number 
density ub by n q = ?i_b/3) retains the ideal gas form 

at Nf r°° 

n q=l J ° 

but with the confinement factor C(T, /j) included. 



In the previous section expressions for the coupling G(T, /i = 0) and the confinement factor C(T, /i = 0) 
are given. These expressions can be generalized to finite chemical potential in a thermodynamically 
self-consistent way using Maxwell relations. Imposing the Maxwell relation between the derivatives of 
the quark number density and the entropy, 



dn 
df 



E 



dn dm 2 



ds dmi 



dm? dT dm 2 d\x 



= and 



dn dC_ 
dCdT 



ds_dC_ 
dC~dfi 



= 0, (10) 



yields a set of first oder quasilinear partial differential equations for the effective coupling constant 
G 2 (T,fi) and the confinement factor G(T, /z): 



dC 2 dC 2 
a T (T, n G 2 ) — + 0/1 (T, u; G 2 )— = b(T, u; G 2 ), (11) 

dC dC 
c T (T^;G 2 )—+c,(T^;G 2 )— = 0. (12) 

The coefficients ar, b, ct, c m depend on T, /i, G 2 but not on C. It can be solved by the method of 
characteristics (see appcndix[5J). The flow of the effective coupling and the confinement factor is elliptic. 
In particular, one finds 

a T (T,fi = 0) = 0, a fj,(T = 0, //) = 0, c r (T, fj, = 0) = 0, c^T = 0, /i) = 0. (13) 

Therefore, the characteristics are perpendicular to both the T and the /i axis. This guarantees that 
specifying the coupling constant and the confinement factor on the T axis sets up a valid initial con- 
dition problem. Plots of the characteristic curves and the confinement factor are shown in figurc^and[5] 



4 Comparison with lattice results 

Simulations of QCD at finite chemical potential are extremely difficult because the fermion determinant 
becomes complex. This prohibits Monte Carlo importance sampling, which interprets the measure as a 
probability factor and thus requires it to be positive. While this problem remains unsolved, there are 
some approaches which circumvent the sign problem and allow lattice simulations for small chemical 
potentials /i <T C . A review comparing these methods in detail can be found in |26) . 
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Figure 1 : Characteristic curves of constant confinement factor C (T, fx) = const, obtained when solving 




Figure 2: The confinement factor C(T,/j,) as a function of the temperature T and the quark chemical 
potential /i. 



G 



4.1 The phase boundary line 



In the case of vanishing chemical potential, universal arguments and lattice simulations suggest a phase 
transition from the hadronic phase to the QGP phase at a critical temperature T c . For QCD with three 
light flavors m u ~ rrid ~ m s ~ 5 MeV this transition is expected to be first order. For two light flavors 
m u ~ rrid ~ 5 MeV and an infinitely large m s there is no phase transition, only a smooth crossover |27|. 
This suggests there is a critical strange mass m c s at which one finds a second order phase transition. 
Lattice calculations indicate that m c s is about half of the physical mass m s . At finite quark chemical 
potential fi and vanishing T a first order phase transition is predicted. For the physical m s this implies 
that there is a first order phase transition for small T and large /x which ends at a critical point (T* , fj,*). 
At this point the phase transition is of second order. For large T and small \i the two phases are separated 
by a crossover. We refer to the line T c {p) that separates the hadronic phase from the QGP phase as the 
"phase boundary line" . In the literature [HI CI IE] this line is also frequently called the "pseudocritical 
line". T c (fj,) has been calculated on the lattice for N f = 4 E| , N f = 2 [7] and N f = 3 flavors 
of quarks up to quark chemical potentials \x of order T c . In the following we focus on the three- flavor 
results where the critical line has been calculated with an accuracy up to terms of order (/i/T) 6 . There, 
a Wilson gauge action and three degenerate flavors of staggered quarks have been employed, with bare 
masses in the range 0.025 < am < 0.04, where a denotes the lattice spacing. The finite volume scaling 
behavior was monitored by using three lattice sizes, 8 3 x 4, I0 3 x 4 and I2 3 x 4. 

In our quasiparticle model, the sudden decrease of the pressure, the energy density, the quark number 
density and the entropy density caused by gluons and quarks getting trapped in glueballs and hadrons 
when T c is approached from above, is parametrized by the confinement factor C(T, //). Consequently, it 
is natural to relate the critical line to the characteristic curve of the confinement factor through T c (/i), 
as long as fi is small and the nature of the quasiparticles does not change qualitatively. 

In order to calculate the confinement factor at finite chemical potential, we need to specify a valid initial 
condition, e.g. C(T, \i = 0). The functional form of C(T, /j, = 0) is set by eq.@. We have employed the 
following set of parameters, as found in ref.|22|: 





Co 


S c 




3 flavors 


1.03 


0.02 


0.2 



We have checked that the form of the phase boundary line in the quasiparticle model depends only 
weakly on the exact choice of parameters and a small difference only shows up for values much larger 
than the range of [i covered by the lattice simulations. The lattice phase boundary line and our result 
is shown in figure|3l The quasiparticle result is within the lattice estimate for [ib < 2.5T C and deviates 
only slightly from the lattice result for larger chemical potentials. 

4.2 Thermodynamical quantities 

There have been lattice calculations of thermodynamical quantities at finite chemical potential for Nf = 
2 + f an d Nf = 2 10 flavors of quarks. In the following we focus on results from JU] where a 
p4-improved staggered action on a 16 3 x 4 lattice was used. There, the N T dependence is known to be 
small, in contrast to standard staggered fermion actions which show substantially larger cut-off effects. 
Estimates of the pressure, the quark number density and associated susceptibilities as functions of the 
quark chemical potential were made via a Taylor series expansion of the thermodynamic grand canonical 
potential Q up to fourth order. 

To calculate thermodynamical quantities within the quasiparticle model, we need to fix the parameters of 
the effective coupling constant and the confinement factor. Our calculations have shown that the results 
are not sensitive to the detailed choice of parameters for the effective coupling G. We have therefore used 
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the parameters from ref. 22 in our calculations. In principle, the parameters of the confinement factor 
can be fixed by comparing our calculations to lattice results at vanishing chemical potential. However, 
in ref.^U] no /x = lattice data is given. Since lattice calculations including quarks give slightly different 
results depending on which action has been used, fitting the parameters in C(T, zx = 0) by comparing 
quasiparticle results to lattice data from a different simulation is not feasible and would lead to large 
differences. Consequently, we directly used the finite /i lattice results for fitting. Good agreement with 
the lattice thermodynamical observables was found for the following sets of parameters: 





Co 


So 


/3 C 


Set A 


1.05 


-0.016 


0.15 


Set B 


1.12 


0.02 


0.2 



While set A reproduces the lattice thermodynamical results slightly better, set B is in better agreement 
with the parameters found in |22j for fi = lattice simulations. 

The temperature dependence of the normalized pressure difference Ap(T, fi) = (p(T, p) —p(T, p = 0))/T 4 
is shown in figure 0] and that of the normalized quark number density n q (T, p)/T 3 in figure Whereas 
the computation of the quark number density from equation is straightforward, a numerical evaluation 
of JBJ is difficult because of the derivatives of the effective masses and the confinement factor in B(T, p) 
(see expressions in appendix lA")l . It turns out that it is simpler to calculate the pressure difference using 
the following relation: 

Ap(T, M ) = ^ £dp,'n q (T,p'). (14) 

The lattice pressure difference is well reproduced even for the largest values of the chemical potential. 
The quark number density is in very good agreement with the lattice data for p/T c — 0.2 and 0.4. For 
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larger values of fi our calculations underestimate the magnitude of the lattice results close to T c , but 
show the same qualitative features. 



5 Summary 

We have presented a quasiparticle description of the QCD EoS at finite temperature and quark chemical 
potential. Our main modification as compared to previous work is the inclusion of finite quark chemical 
potential in a thermodynamically consistent way. We have first reviewed our improved quasiparticle 
model which schematically includes confinement. We have shown how Maxwell relations can be used to 
construct the effective coupling G(T, fx) and the confinement factor C(T, fx) at finite chemical potential. 
We then used this model to calculate the phase boundary line T c (fx) and the normalized pressure differ- 
ence Ap(T,(x) = (p(T,fx) — p(T,fx — 0))/T 4 and the normalized quark number density n q (T, lm)/T 3 . We 
compared our results to recent lattice calculations and found remarkably good agreement even for large 
quark chemical potentials fi ~ T c . 

A Calculation of B(T,fj) 

The "background field" quantity B(T,fi) appearing in eq.ljSJ) can be obtained from the Gibbs-Duhem 
relation 

e+p = Ts + fxn = T—+fx—. (15) 

The left hand side reads: 

N r Nf r°° , /4fc 2 + 3m 2 \ 

?-+P=i^J o dk[f+ + f D ]C(T,fx)k 2 ( ^ q \ . (16) 

To evaluate the right-hand side, derivatives of f^(E^) with respect to T and fi are rewritten as derivatives 
with respect to k. After an integration by parts, the first term on the right-hand side reads 

Bp N c Nf r„ f+ (dCtf 3k 2 dm 2 h* \ 

= T ^T j Q dkf D\-Qf^- k ~ C ( T ^) ? ^ fyr +C(T,fx)(E k + fx)+C(T,fx) 



dT 3tt 2 J JU \dTE k v 2E k dT v f/ v TE k I 

+ T M± rdkrJ^^CiT^^^ + CiT^m-f^CiT^)^ 
3tt 2 J JD ydTE k y 2E k dT y ' WV k w V TE k J 

- (17) 
dT ' { '> 

and the second term is given by 

Substituting (|16|) , (|17|) and l|18|l in the Gibbs-Duhem relation yields a partial differential equation of the 
type 

df(x,y) df{x,y) 
x—^ + y—^ r =l(x,y). (19) 
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It has the general solution 

f(x,y)= P dtl(t,h)+n(^). (20) 

J x VI/ 

Here, TL{y/x) is a solution of the homogeneous equation. Returning to our case, W(/x/T) becomes an 
arbitrary function of the ratio n/T to be fixed by boundary conditions. For // — > 0, 7i(/x/T) does not 
depend on T anymore and therefore has to be identified with an integration constant Bq. Provided 
that Tt(/ji/T) is a continuous function it must be close to Bo for small n/T. The first term in a Taylor 
expansion of H.(fi/T) vanishes and the series starts only at order (/j,/T) 2 . Therefore we identify 7i(n/T) 
with the constant Bq for all /i under consideration. Assembling all pieces, the final result reads 

B(T,fi) = B 1 (T, fJ ,) + B 2 (T, fJ ,)+B , 

where the explicit r-dependence in C(r, /i/Tt), m q (r, /i/T t) and E%(t, ji/T r) has been suppressed for 
the sake of lucidity. 

B Method of characteristics 

Equations (|llfl and i|12fl are a set of coupled quasilinear first order partial differential equations for the 
effective coupling constant G 2 (T, /i) and the confinement factor C(T, (i). Equation does not depend 
on C(T, /x). Thus we can first solve this equation for G 2 (T, /x) and insert the result in equation 112(1. 



The usual method found in textbooks is to reduce a quasilinear partial differential equation of the form 

— +a,(T, A ,;X) — 
to a system of coupled ordinary differential equations, 



a T (T,fi;X)— + a^T,ti;X)— = c(T>;X) (22) 



dT( S ) d M ( S ) dX( s ) 

^ = flT ' "dT^' (23) 
This determines the characteristic curves T(s), /i(s), and the evolution of X along such a curve, given an 
initial value. However, this method is not well suited for numerical use which is necessary for non-trivial 
ax, a M and c. Rewriting equation (1221 as 



dX dX d/i \ dX dX ( a 



c 



a T{w-^) +a ^= c ^^ [t dT -v=- T dT - dx ^ (24) 

we find the equation a^dT — ard/i = for the characteristics and cdT — ardX = for the evolution of 
X. These equations can easily be solved numerically. 
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